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Any viable stellarator reactor will need to be nearly omnigenous, meaning the radial 
guiding-center drift velocity averages to zero over time for all particles. While omni- 
genity is easier to achieve than quasisymmetry, we show here that several properties 
of quasisymmetric plasmas also apply directly or with only minor modification to 
the larger class of omnigenous plasmas. For example, concise expressions exist for 
the flow and current, closely resembling those for a tokamak, and these expressions 
are explicit in that no magnetic differential equations remain. A helicity (M, N) 
can be defined for any omnigenous field, based on the topology by which B contours 
close on a flux surface, generalizing the helicity associated with quasisymmetric fields. 
For generalized quasi-poloidal symmetry (M = 0), the bootstrap current vanishes, 
which may yield desirable equilibrium and stability properties. A concise expression 
is derived for the radial electric field in any omnigenous plasma that is not quasisym- 
metric. The fact that tokamak-like analytical calculations are possible in omnigenous 
plasmas despite their fully-3D magnetic spectrum makes these configurations useful 
for gaining insight and benchmarking codes. A construction is given to produce 
omnigenous B(9, Q patterns with stellarator symmetry. 
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I. INTRODUCTION 



Nonaxisymmetric toroidal magnetic fields can provide intrinsically steady-state, disruption- 
free plasma confinement. However, unlike axisymmetric fields, nonaxisymmetric fields do 
not generally confine all trapped particle orbits. This shortcoming has a modest deleterious 
effect on energy confinement, which still follows scalings comparable to tokamak confinement 
due to turbulent transport^. However, unconfined orbits will likely pose a serious problem 
in a reactor, where unconfined energetic alpha particles will collide with the first wall before 
thermalizing, causing unacceptable damage to the plasma-facing components. Consequently, 
one criterion of stellarator optimization in recent years has been the quality of collisionless 
particle confinement. The ideal limit in which all collisionless trajectories are confined is 
known as omnigenity (sometimes spelled omnigeneity) . Omnigenity can be defined more 
precisely as the condition that the time average of v m -Vip along each field line vanishes for all 
values of magnetic moment /i = vj_/(2B). Here, v m = v^B x K,/(£lB)+vj_(B x VB)/{2VtB 2 ) 
is the magnetic drift velocity, n = b- Vb, b = B/B, B = \B\,Q = ZeB/(mc) is the gyrofre- 
quency, Z is the species charge in units of the proton charge e, m is the mass, c is the speed 
of light, and 2irip is the toroidal flux. As shown in Appendix an equivalent definition 
of omnigenity is that the longitudinal adiabatic invariant J = §vud£ is a constant on a 
flux surface. Other equivalent definitions are the absence of a "1/V" regime of confinement, 
where v is the collision frequency, or an effective helical ripple of zero. 

One method of obtaining omnigenity in a nonaxisymmetric toroidal system is quasisym- 
metry, the design principle behind the HSX- and NCSX- experiments. Quasisymmetry is 
usually defined^ as the condition that the field magnitude B varies on a flux surface only 
through a fixed linear combination of the poloidal and toroidal Boozer angles: 

B = B(ip,M6-N() (1) 

for integers M and N. Quasisymmetry can also be defined as B = B(ip, M9*—N£*) for other 
coordinates (#*, £*) such as Hamada angles (as proved in Appendix [B]), or by the coordinate- 
free conditions^ B ■ V [(B x Vip ■ VB)/(B ■ VB)] = or-& VBxV^p- V(B ■ VB) = . The 
Lagrangian L for guiding-center drift motion, when expressed in Boozer coordinates^, only 
depends on B and not the vector components of B. Consequently, an ignorable coordinate in 
B gives rise to a conservation law. The conserved quantity^, resembling canonical angular 
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momentum, ensures each particle can only drift a distance on the order of a gyroradius away 
from a given flux surface, implying omnigenity. 

While every quasisymmetric field is omnigenous, not every omnigenous field is quasisym- 
metric. A proof by Cary and Shasharina at first appears to show the opposite conclusion, 
that a field which is both infinitely different iable (analytic) and perfectly omnigenous must 
be quasisymmetrio^ 1 ^. However, the same authors point out that the proof is quite fragile: 
an analytic field that deviates from omnigenity only slightly may depart from quasisymme- 
try by a great amount. (We will construct examples of such fields in Section [V]) Thus, in 
practice, omnigenity does not imply quasisymmetry. Consequently, to achieve good colli- 
sionless particle confinement, there is no need to strive for the strict condition of quasisym- 
metry when the weaker condition of omnigenity is sufficient. While it is possible to find 
three-dimensional MHD equilibria that are reasonably close to quasisymmetry^, lifting the 
demand of quasisymmetry widens the parameter space, allowing stellarators to be better 
optimized for other criteria. 

One may still aim to achieve quasisymmetry for a different reason: In non-quasisymmetric 
plasmas, the parallel flow is determined by leading-order neoclassical processes, while in 
quasisymmetric plasmas, the parallel flow is determined by turbulence^. As flows and flow 
shear may affect MHD modes and microinstabilities, quasisymmetric plasmas may have 
unique stability and turbulence properties, though more research is required to explore this 
notion. 

Quasisymmetric fields provide an important and useful ideal limit for understanding non- 
axisymmetric plasmas. In quasisymmetric fields, concise analytic expressions can be derived 
for the neoclassical distribution function, radial fluxes, and parallel flows and current&ii. 
These formulae are nearly identical to the corresponding formulae for a tokamak. The for- 
mulae are explicit, in that they do not involve solutions of partial differential equations. In 
contrast, the corresponding formulae for a general stellarator must be expressed in terms of 
the solutions of partial differential equations involving the field strength^— . One example 
is the system of equations (I571)-( l60|) for the bootstrap current. 

It is the purpose of this paper to demonstrate that perfect omnigenity is another useful 
ideal limit. Omnigenity, while less restrictive than quasisymmetry, is still a strong enough 
condition to place powerful constraints on B(6,()- For example^, in the neighborhood of 
each minimum and maximum of B along a field line, each adjacent field line segment on the 
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flux surface has an extremum at the same value of B, as we will show in Section |TT] using 
novel arguments. We will also show that field lines can never be parallel to the contours of B 
on each flux surface, and these contours must topologically link the flux surface toroidally, 
poloidally, or both. By defining integers M and iV as the number of times each contour 
links the plasma toroidally and poloidally respectively, the helical "mode numbers" M and 
N in (pQ) can be generalized to any omnigenous plasma. Using this generalization, we 
will show in Section III II that formulae for the bootstrap current and flux-surface-averaged 
parallel flow in a quasisymmetric plasma in fact apply to the larger class of all omnigenous 
plasmas. The formulae for the Pfirsch-Schliiter flow and current in a quasisymmetric plasma 
also apply to all omnigenous plasmas if a new term is added, consisting of an integral of 
derivatives of B. The same is true of the low-collisionality distribution functions. Some of 
these neoclassical properties were discussed previously in Refs. 
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for the specific 



case of generalized poloidal symmetry (M = 0). Here we give more general calculations that 



apply to omnigenity of any helicity. In Refs. [L6| andll7| it was pointed out that for an M = 



omnigenous plasma, the bootstrap current vanishes unless there is inductive or RF-driven 
current. This result will be recovered from the analysis of general- helicity omnigenous fields 
here. 

Another interesting feature of omnigenous magnetic fields that we will demonstrate con- 
cerns the radial electric field E r . In a general stellarator, the radial electron and ion particle 
fluxes exhibit different dependencies on E r . Consequently, it is possible to determine the 
electric field using the condition that the net radial current must vanish, i.e. quasineutrality 
or ambipolarity. However, this procedure fails in tokamaks or in quasisymmetric stellara- 
tors, which possess the property of "intrinsic ambipolarity." This is the property that to 
leading order in the expansion of gyroradius to system size, the net radial current vanishes 
regardless of the radial electric field, and so the electric field is determined by higher-order 
processes that are more difficult to calculate. An omnigenous plasma that is far from qua- 
sisymmetry will not be intrinsically ambipolar, and so it is still possible to solve for E r using 
ambipolarity. In Section HVl we will derive E r for any omnigenous, non-quasisymmetric field. 
The result differs from previously known expressions for the electric field in non-omnigenous 
stellarators. Our calculation will also show from a new perspective that E r indeed becomes 
undetermined in the limit of quasisymmetry. 

Finally, in section [V] we will derive some further geometric properties of B(6,Q f° r an 



omnigenous flux surface. It will be shown how to generate a family of omnigenous flux 
surfaces that are consistent with an additional symmetry usually possessed by stellarator 
experiments. 

II. MAGNETIC FIELD PROPERTIES 
A. Extrema of B on each flux surface 

Let us now begin to analyze the geometric properties of omnigenous fields in detail. Let 
B(r) and B(r) denote the nearest minimum and maximum of B found by moving forward 
and backward along a field line from the starting position r. In a general stellarator, B and 
B will take on a continua of values (or several continua) on each flux surface. However, in 
an omnigenous field we will now show that B and B may only take on discrete values on 
each flux surface, and in the simplest case, each has only a single value on each flux surface. 
In other words, in the neighborhood of an extremum of B along a field line, each nearby 
field line segment on the same flux surface has an extremum at the same value of B. 

Let us first prove this property for the minimum B. We begin by recalling the con- 
travariant and covariant expressions^ for B in terms of the poloidal Boozer angle 9 and the 
toroidal Boozer angle (: 

B = Vip x V0 + *VC x V^, (2) 

B = 0(i>, 9, C) W> + I(i>)V9 + Gty)VC. (3) 

Here, * is the rotational transform, is 2/c times the toroidal current inside the flux 

surface ip, and G(tp) is 2/c times the poloidal current outside the flux surface ip. Now let 
r o — {9qi Co) denote a point at which B is minimized with respect to movement along a field 
line. At this point, B ■ VB must vanish, so 

= (B ■ V0) (dB/dd) + (B ■ VC) (dB/dC) ■ (4) 

If the field is omnigenous, deeply trapped particles at ro must have no radial drift. As 
v m ■ Vtp oc B x VB ■ Vijj, then 

= B x VB ■ Vip = —G(B ■ V()(dB/d9) + I(B ■ V()(dB/d(). (5) 

Equations (@| and flS} are a system of two linearly independent equations for dB/d9 and 
dB/d( at r , so both of these derivatives must be zero. The vanishing of these derivatives 
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implies B(vq) is either an isolated local minimum, a saddle point, or one point along a 
"valley" of constant B. 

The first two of these three possibilities can be excluded as follows. We first Taylor-expand 

B*B + ^(9- 9 ) 2 + y(9 - fl„)(C - Co) + f (C - Co) 2 , (6) 

where x = [d 2 B/d9 2 ] , y = [d 2 B/d9d(] , z = [d 2 B/d( 2 } , and the zero subscripts indicate 
quantities evaluated at Tq. Considering the variation of B along the nearby field line 9 = 
$o + *(C — Co) + $ f°r some small 5, by eliminating either 9 or C in (JS]) and completing the 
square, it is evident that B is minimized along the field line at the point r*i = Ci), where 
9i = 9 + (-by + z)6/A, Ci = Co ~ ( iX + y)8/A, and A = t 2 x + 2hy + z. Plugging the definitions 
of 9\ and Ci into (jSJ) gives 

B^B + ^ + (9-9 1 )^-(C-C 1 )^ + ^9-9 1 ) 2 + y(9-9 1 )(C-C 1 ) + l(C-^ 2 (7) 

where k = xz — y 2 . Since B is minimized along the shifted field line at 7*1, the radial drift 
must vanish there, for the same reason it had to vanish at t*o, so dB/d9 = and dB/d( = 
at r*i. Consequently the third and fourth right-hand-side terms in (J7J) (those linear in {9 — 9\) 
and (C — Ci)) niust vanish, and so k must be zero. Thus, the second term on the right hand 
side of (JTj) vanishes, and so the minimum of B along the shifted field line is Bq, the same as 
the minimum on the original field line. This proves the desired result for B. This fact, that 
deeply trapped particles are perfectly confined when all B on a flux surface are the same, 
was first observed in Ref. 20. 



Now we make the analogous argument for B. Let r now represent a point at which B is 
maximized along a field line. The bounce time diverges for barely trapped particles, because 
they spend an infinitely long time near the turning points. Thus, if there is an outward radial 
drift at r , marginally trapped particles there will have an arbitrarily large radial excursion. 
Even if these particles would in principle make a large inward radial step at the opposite 
bounce point, they would drift out of the machine before having time to get to the other 
bounce point, so these particles effectively would have a nonzero time-averaged radial drift. 
We choose to include in the definition of omnigenity the condition that marginally trapped 
particles cannot have radial steps of unbounded size in this manner. Now suppose at vq 
there were an inward radial drift. By omnigenity, marginally trapped particles must make 
an arbitrarily large outward excursion elsewhere in the trajectory to balance the arbitrarily 
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FIG. 1. (Color online) Omnigenous fields have an intermediate level of complexity between qua- 
sisymmetric and general nonaxisymmetric fields. 

large inward excursion in the neighborhood of r , so this case too is unacceptable. Therefore, 
v m ■ Vip must vanish at t*o, implying (jSJ). The rest of the argument for the constancy of B 
then applies, and so B must be constant as well for each field line in a neighborhood of r*o 
on the flux surface. 

Due to these constraints on the extrema of B, omnigenous fields represent an intermedi- 
ate level of complexity between quasisymmetric fields and general nonaxisymmetric fields. 
Figure [1] illustrates this point. For axisymmetric and quasisymmetric fields, the B wells in 
which particle are trapped all have identical shape. At the opposite extreme of a general 
three-dimensional field, different field line segments have different maxima and minima of 
B. Omnigenous fields represent a middle ground. Maxima of B occur repeatedly at the 
same values of B, and the same is true of the minima. In fact, we will prove in Section Til CI 
that the maxima are equally spaced (in 6, in (, and in distance along the field), as indicated 
by the horizontal red arrows in Figured] Yet, unlike in quasisymmetric fields, the shape of 
the B wells is different at different points along a field line. 

B. Generalized helicity 

In an omnigenous toroidal field, the constant-i? contours on a flux surface must all encircle 
the plasma toroidally, poloidally, or both. In other words, each contour must topologically 
link the flux surface: a contour cannot be continuously deformable (homotopic) to a point 



FIG. 2. (Color online) Contours of B on a flux surface of a non-omnigenous stellarator. Contours 
typically exist that do not topologically link the flux surface, such as the bold contour here. The 
bold black straight line is a field line. 

without leaving the flux surface. If any constant-5 contour did not link the flux surface in 
this manner, then the contour would enclose a point maximum or minimum of B within 
the surface, and we just proved that such point extrema cannot occur in an omnigenous 
field. As an illustration, the bold curve in Figure |2] is a B contour that does not link the 
plasma. This contour encloses a point minimum or maximum of B(9, () at P, and such 
points cannot exist in an omnigenous field. For contrast, Figure [3] depicts an omnigenous 
field, one in which the B contours encircle the plasma both poloidally and toroidally. Figure 
H] shows an omnigenous field with toroidally closed B contours. Neither field has any point 
extrema of B. 

We can also prove that all B contours must encircle the plasma using the following 
alternative argument. If a constant-I? contour does not link the flux surface, then there will 
be points at which the contour is tangent to the field, such as point T in Figure El (Recall 
that field lines are straight in the (6, Q plane for Boozer coordinates.) At such a point of 
tengency, B ■ VB = 0, but the derivative of B in any other direction on the flux surface 
is nonzero. Thus, (jlj) is satisfied while ([S]) is not, violating the condition of omnigenity. 
Physically, if B is a minimum along the field line at T, deeply trapped particles at T will 
have nonzero average radial drift, while if B is a maximum along the field line at T, barely 
trapped particles will make an infinite radial step when they bounce at T. Neither type of 
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FIG. 3. (Color online) The field magnitude B / (B 2 ) 1 / 2 for two of the four periods of an omnigenous 
field with M = 1, N = 4, and * = 1.05. The straight dashed line is the maximum and the curved 
dotted curve is the minimum. Field lines are parallel to the arrow. One branch is shaded and the 
other is unshaded. 




FIG. 4. (Color online) The field magnitude B/{B 2 ) 1 ^ 2 for two of the three periods of an omnigenous 
field with M = 1, N = (generalized quasi-axisymmetry) , and * = 1.62. Markings are as defined 
in Figure [3l The arrows point along B and connect two contours of equal B. The conditions 
dAg/dx = and dA^/dx = imply the arrows must have the same length. 
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radial motion is allowed in an omnigenous field, so the B contours of each flux surface in an 
omnigenous field can never be tangent to field lines. This implies the contours must link the 
plasma. It can be seen in Figures [3] and @] that the B contours are indeed nowhere tangent 
to the field lines. 

We can now define integers M and N as follows: each B contour closes on itself after 
traversing the torus M times toroidally and N times poloidally, that is, after 9 increases 
by 2irN and £ increases by 2ttM. This convention may seem backward at first, but it is 
consistent with the M and N defined earlier for a quasisymmetric field: B = B(M6 — N(). 
By defining M and N in terms of the topology of the B contours as we have done, the 
helicity associated with quasisymmetric fields is generalized to any omnigenous field. Notice 
that this helicity is completely independent of the rotational transform 

It should be noted that definition of the term "quasi-iso dynamic" given by some authors^ 
is equivalent to the condition "omnigenous with M = 0." However, other authors^- define 
"quasi-isodynamic" differently, as the case in which only particles with a particular value of 
normalized magnetic moment A = v\/{v 2 B) are omnigenous, not all particles. 

Next, it will be convenient to define a field line label x by 

X=(0-<)/(iV-*M). (8) 

This definition is convenient because if a constant-l? curve is followed until it closes on 
itself, then x wn l increase by 2tt. For much of the analysis that follows we will use 
(ip, x, B) coordinates. In these coordinates the magnetic field has a contravariant form 
B = (N — *M) Vijj x Vx and a covariant form 

B = B^Vip + B B VB + B X V X - (9) 

The Jacobian is x Vx • V5 = (B ■ VB)/ (TV — tM). We now derive several properties 
of the covariant coefficients. The inner product of (Q with B gives Bb = B 2 / B ■ VB. The 
inner product of (Q with x VB gives 

B X = -(N-*M) . (10) 

As discussed in Section Hi Al the numerator of ( TIOl) vanishes whenever the denominator does, 
leaving B x nonsingular. Indeed, since omnigenity precludes B contours from being tangent 
to field lines, then dr/dx is never singular, and so B x = B-dr/dx cannot be singular. Here 
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FIG. 5. (Color online) Equation (jlip is derived by applying Ampere's Law to a B contour, such 
as the red curve here. The black arrows illustrate the toroidal and poloidal currents, which are the 
currents through the red and blue translucent surfaces respectively. 

and throughout this paper, d/d\ is performed at fixed B, 3/3B is performed at fixed x> 
3/39 is performed at fixed £, and 3/3( is performed at fixed 9, unless denoted explicitly 
with a subscript. 

Now consider a path on a flux surface that follows a constant-I? curve until it closes 
on itself. As described above, x increases by 2tt along this curve. From Ampere's Law, 
(4ir/c)ih = JB ■ dr — Jq W cZ%B ■ 3r/3x = Jo" B x dx where ih is the current linked by the 
loop, and the integrals are performed at constant ip and B. The integration path links 
the torus N times poloidally and M times poloidally, so the amount of current linked by 
this helical curve is ih = Ni t + Mi p where it is the toroidal current inside the flux surface 
and i p is the poloidal current outside the flux surface. These currents are related to the 
coefficients of the Boozer covariant representation by I = 2i t /c and G = 2i p /c. Therefore 
i h = c (MG + NI) /2 and J 27r B x dx = 2ir (MG + NI). This application of Ampere's Law is 
illustrated in Figure [5] (for M — 1 and N — 4.) Any single- valued function of position must 
be periodic in x (with period 2tt) if B is held fixed. In particular, B x must be periodic in 
this way, so we can write B x = MG + NI + 3h/3x for some single- valued h. Hence, recalling 
f fTUl) . we obtain the useful formula 



In the quasisymmetric limit B(8,Q = B (M8 — N(), then the left-hand side of this 




(11) 
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expression can be computed directly, and the result is the same but without the h term. 
Thus, quasisymmetry corresponds to the dhjd% — limit. 



C. Relation between branches 

The continuous coordinates (if), x, B) do not uniquely determine a point in space, both 
because there may be multiple toroidal segments, but also because within each segment there 
are two points at given B on either side of B, the minimum of B on the flux surface. This 
discrete degree of freedom, called the "branch," will be denoted by 7 = ±1. The shaded 
and unshaded regions of Figures |3] and H] illustrate the two branches. The variations of B in 
the two branches are related due to the condition of omnigenity. As shown in Appendix [Aj 

^E^— = 0. (12) 
d x V b • vs 



Theorem" in Ref. 
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Using (Tl2l) . we can prove several facts about pairs of points on a same field line that 
share the same B but lie on opposite sides of B. Let be the distance between these 
points, measured along the field line. In a general stellarator, Ai will depend on the field 
line label x (in addition to B and ip). But since A e = E 7 7 ■ VB)'Y 1 dB' (i.e. b ■ VB 
is evaluated at B' rather than B), then ( fl2|) implies dAe/dx = in an omnigenous field. 
Thus, for any given B, every such pair of points has the same separation Af . This result is 
illustrated for the case of B, the maximum of B on the flux surface, by the horizontal red 
arrows in Figure [TJ A similar result holds for the difference in ( between pairs of points. 
Let A^ = J2-yl fi ' (d(/dB)'dB' be the difference in ( between a pair of points as described 
above. Notice d(/dB = (B ■ \7()/B ■ \7B. Multiplying the covariant and contravariant 
Boozer representations (E])-©, 

B 2 /(G + hi) = Vif) x V6> ■ VC = B • VC = * _1 B ■ V0. (13) 

Therefore, 



1 f B b dB 

^ = gT^? 7 7b 7*TW (14) 



which must be independent of x due to (fl2"j) . A similar proof shows the separation in 9 
between the points is also independent of field line. The results dA 9 /dx = and dA^/dx = 
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are illustrated by the three arrows in Figure HI These arrows point along B, all joining two 
contours of the same B on opposite sides of B. As dAg/dx = and dA^/dx = 0, these 
arrows must all have the same length. 

It follows that the contours of B = B (where B is again the maximum B on the flux 
surface) must in fact be straight lines in the [6,() plane. The basis of the argument is that 
as Af(jB) and Ag(B) are constants on a flux surface, then when the B contour is translated 
by Aq(B) and Ag(B), it must lie on top of itself. In other words, the B contours must be 
symmetric under a translation along B, shown by the arrow in Figure [31 Except for the 
uninteresting case in which * is a special low-order rational number, the B contour cannot 
possibly have this symmetry unless it is straight. The rotational transform * can be assumed 
to be irrational, since by continuity, B on any rational surface can differ only infinitesimally 
from Bona nearby irrational surface. To begin the rigorous proof, first consider the M = 
case (poloidally closed B contours), and suppose the stellarator has N p identical toroidal 
segments with one B contour per segment. Let one B contour be given by £ = y{6). If 
we shift this contour by Ag(B) and A$(B), it must lie on top of the next B contour, the 
one given by C = y(0) + (2tt/N p ). Therefore, y(0 - *A C ) + A c = y(6) + (2n/N p ). Then 
expanding y as the Fourier series y(9) = J2 n 34e m6 \ we can write 

^»e w (l-e^<)=A{-p (15) 

The n = component of this equation implies A^ = 2n/N p . It follows that the exponent 
—imA^ will never be an integer multiple of 2ni if * is irrational. Therefore, the quantity in 
parentheses in f JT5|) can never be zero for n^O. Every n^O component of ffl5l) consequently 
implies y n = 0, so y(9) must be constant, meaning the B contours are straight. To apply 
the proof to fields in which both M and iV are nonzero, all that is needed is to redefine 
y(0) as £ — M8/N along the B contour so 3^ remains periodic in 0. The proof can also be 
adapted to the iV = case by switching the roles of 9 and £■ 

Finally, it is important to consider both branches when forming the flux surface average 
in the (if), x, B) coordinate system. For any quantity Q, this average is given by 

(Q) = Y^C d *L* dB B%B ^ 
where V = E 7 7 C d X f§ dB/B ■ VB. 
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D. Departure from quasisymmetry 



For the remainder of this section we develop several properties related to dh/dx, a 
quantity that represents the departure from quasisymmetry. These properties will gen- 
eralize results discussed in Ref. (— ). First, plugging into the MHD equilibrium relation 
= Vtf} ■ V x B, we find dBs/dx = dB x /dB. Plugging in our earlier expressions for Bb 
and B x , then 

(d/d X )(B 2 /B ■ VB) = d 2 h/dBd X - (17) 

Applying X) 7 7 an d recalling f|T2|) . then d 2 h/dBdx = 0. Now integrate this expression 
in B from B to B. As dh/dx is continuous everywhere, it is continuous at B, and so it must 
be branch-independent at B. Therefore the contribution from the integration boundary at 
B vanishes. Consequently, 

^1 dh/dx = 0. (18) 

In other words, dh/dx is branch-independent everywhere. 

Any quantity that is continuous and branch-independent must be a constant along the 
curve B = B. This result applies in particular to dh/dx, f° r which the constant must be 
zero, or else J 27r dxdh/dx would be nonzero. Thus, dh/dx must vanish along B = B. 

Now we derive expressions to relate the new quantity dh/dx to more familiar Boozer 
coordinates. Using f fT3|) . 

* iG+lI) (°4 +i °4)-\ (19) 



B ■ VB v ' V dC, d9 
In addition, using (d£/d9) x = a -1 gives 

(dB\ dB ldB 

{m) x = w + :tk' (20) 

where subscripts on partial derivatives indicate quantities held fixed. Combining this result 
with (US]), 

* e + (G+t/) |< (21) 



B ■ VB * dB v ' dB 
Next, we form BxVf VB= (V^ X ■ VC) \G{dB/d6) - I(dB/d()}. Combining this 
result with ( IT9l) gives 

BxVf VB _ ( dB dB\ (dB dB\ l 

b-vb ~ r^" J acJ {~dc + "^e) ' { ) 
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FIG. 6. (Color online) The departure from quasisymmetry, dh/dx, normalized by G + */, and 
calculated for the field of figure [3] using (|23j). Dashed contours are negative. 



Substituting for the left-hand side using (ITT]) , after some manipulation we obtain 



dh 

dx 



-{G + */) 



(23) 



This expression allows explicit calculation of dh/dx f° r a given B(9,Q. As an example, 
Figure |6] shows dh / dx calculated using this formula for the example field of Figure El 

Lastly, we derive some additional formulae for dh/dx that will be used later. From 
dx/dQ = 1/(N - iM), we find dB/dO = (dB/d6) x + (N — hM)~ 1 {dB /d X )e- Plugging this 
result into (123]) gives 



dh 
dx 



(G + */) 



N ■ 



{dB/d X ) t 
(dB/d6) 



{G + a) (de ^ 

* \dx ) 



(24) 



where the second equality follows from = (dB/d8) B = (dB/d8) x + (dx/ d6) B {dB / dx)e- 
Applying d/dx to (jHJ), we find d6 /dx — N = i(d(/dx) — *M. When inserted into ( 1241) . this 



gives 



9h T , fd( ' 



(25) 



Comparing either fflty or ( 125]) with ( 121]) . it is evident that (1T7J) is satisfied. 

Expressions (123]) (125]) will be used in section ITlI Al to calculate the Pfirsch-Schliiter flow 
and current. 
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III. PARALLEL FLOW AND CURRENT 



We now move on to analyzing the physical properties of omnigenous plasmas. First we 
will derive general properties of the flow and current for any collisionality. Then we will 
calculate the flow and current explicitly for the long-mean-free-path regime. 

A. Form of the flows and current 

Consider the density n and flow velocity V of a single species. We take n to be a flux 
function to leading order, and take the perpendicular flow to be given by the sum of E x B and 
diamagnetic flows. Then the flow must satisfy = V-V = B-V(V\\/B)+u~B x Vip-V(l/B 2 ), 
where oj = c(d§ /dtp) + c(dp/dip)(Zen)~ 1 , $ = ($) is the flux surface average of the 
electrostatic potential <3>, p is the species pressure and a flux function to lowest order, and 
we have used the fact that V ■ (B x V?/>) = for MHD equilibrium. Now define U to be the 
single-valued and continuous solution of 

B-V((7/B)=BxVf V(l/5 2 ) (26) 

with the integration constant chosen so that (UB) = 0. The solubility condition of ( 12 6 p is 
satisfied for any MHD equilibrium. Then = B ■ V[(Vj| + uU)/B\, so V\\ + ull = A (ip) B 
for some unknown A, which can be eliminated by multiplying the last equation by B and 
flux surface averaging. Thus 

To solve for U, we use (JTTj) , ([25)1, and B • V = (B ■ VB)d/dB to obtain 

±(°) = * 1 ( MG + NI + ¥). (28) 
dB\Bj B 3 (N-iM)\ d X J 

Integrating, 

U_ _ _ 1 (MG + NI) _ 2 rB dB' dh' 
B~ ~fP (N--tM) ~ (N- *M) Jb (B'f d% 
where Y(ip) is an integration constant. The limit of integration has been chosen to ensure Y 
is continuous, which can be seen as follows. The parallel flow Vn has a term proportional to 
U (in f)27p ). so U must be continuous everywhere. At B = B, then U must be independent 
of 7, and so from ( )29|) . Y then must be independent of 7 (since the other terms are all 
7- independent.) Next, consider that aside from perhaps Y, all the other terms in ( 129]) are 
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continuous across the curve B = B at fixed 9 (or, in the N = case, at fixed £•) Therefore 
Y must have this same property. Therefore Y must be independent of %■ 

The unknown Y can be determined by multiplying (129 j) by -B 2 and flux surface averaging. 
Defining 

^ d5' 



W(ifj lX ,B) = 2B 2 f 
Je 



b {B'fdx 



then 



B (N — *M) 



(MG + iV/) ^1 



5 s 



(30) 



(31) 



Since we constructed U to satisfy (UB) = 0, it should be the case that (WB) = 0. 
Indeed, it can be verified that (WB a ) = for any a as follows. First, using (|T2|) and ([18]) . 
we find 



r 277 1 



E 



7 



2w d x —] 







(32) 



VBd X 

where the last equality follows because the last term in parentheses is zero. Examining the 
form of the flux surface average in (fl6l) . it can be seen that (WB a ) is proportional to the 
leftmost expression in (|32|) . so it vanishes. 
Using fl3Tl) . the parallel flow is 



In 



where 



Y 



PS 



(V\\B)B 
(B 2 ) 



1 dp 



PS 



(33) 



B- 



(B 2 ), 



(MG + NI) + W 



B(N -iM) \ dip Zen dip j 
Summing over species a to form j\\ = T^ a Z a en a V\\ a , and using quasi-neutrality, 

, U\\B)B 



(34) 



PS 



(B 2 ) 



where 



PS 



c ( dps 
B (N — iM) \ dip 



B 2 \ 
W)J 



(MG + NI) + W 



(35) 



(36) 



and ps is now the total pressure J2 a Pa- Notice that ( 133]) - (136]) are valid for any collisionality 
regime. 

In a general stellarator, the Pfirsch-Schliiter flow and current can only be defined implic- 
itly, in terms of the solution of the magnetic differential equation ( 126]) . In omnigenous fields, 
in contrast, the flow and current may be written explicitly, in terms of the integral of the 
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FIG. 7. (Color online) The quantity W that arises in the Pfirsch-Schliiter flow and current, nor- 
malized by (G + */), and computed for the magnetic field in Figure El 



field strength W. We now describe two approaches to efficiently calculate W for a given 

B(e,c). 

In the first approach, the integration variable in ( 130]) is changed to 9 using ( 120|) . and then 
(123]) is applied, giving W = 2B 2 (G + tl)w/t , where 

w = / o M— — + N—— = * / — ^ M— - + iV— - . (37) 

Here, 9 and ( represent the values of 9 and £ associated with B = B. Notice the integrals 
in f )37|) are evaluated along constant-x paths (field lines). 

In an alternative approach, (|24|) or ( 125]) is substituted into ( 130]) to obtain 

Numerical root-finding can be used to compute ((x, B) or 9(x, B) from a given B(9, £), and 
the result used to compute either of the integrals in (138]) . Figure [7] shows the W computed 
by either of these methods for the magnetic field of Figure [3j 
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B. Long-mean-free-path regime 



Explicit formulae for (V\\B) and (j\\B) can be derived in the long-mean-free-path regime. 
We begin with the drift-kinetic equation 

v„b ■ Vh + v m ■ + ■ V$x = C {fx} , (39) 



where C is the linearized collision operator and $1 = $ — $o- m f|39|) and hereafter, the 
independent velocity-space variables are A = vj_/(Bv 2 ) and the leading-order total energy 
mv 2 /2 + Ze$ - 

It is useful to next define A by the relation 

v m ■ Vtp = v\\b ■ VA. (40) 

The radial magnetic drifts must have such a form because, from the definition of omnigenity, 
v m ■ V"0 vanishes upon a transit or bounce average. Observing B x k- Vip = b x V-B • Vip, 
then v m ■ Vip = —(B x VB ■ V^)){v\\/ B)d / dB(v\\/Vt). Then ([TTJ and (00]) imply 

|* 1 ( MG + NI + d ±\pX (41) 

dB (N — tM) \ dx) OBVt K 1 

Integrating in B, we find 



l /*. JD/ 0/i' a «f. 



where 

5 = / cLB'^-^^, (43) 

N-tMJB dxdB'tt'' V ; 

and = I? if A < 1/B and -B x = 1/ A otherwise. The limit of integration is chosen this way 
for passing particles (A < 1/B) so that A is continuous in position space across the curve 
B = B, and the limit of integration for trapped particles (A > 1/B) is chosen so that A is 
continuous in velocity space at A = 1/B. 

Thus, the kinetic equation may be written v\\b ■ Vg = C{f\}, where 

9 = / 1 + A§ + ^/.. (44) 

For low collisionality, we expand g = g^ + + . . . and take the leading-order kinetic 
equation to be v\\b ■ Vg^ = 0, so dg^ /dB = 0. For passing particles, continuity of g^ at 
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B requires to be constant on a flux surface, whereas for trapped particles, may still 
depend on x- The next-order kinetic equation is 

«l|6.V^ = C{^-Afl/ o /a0}. (45) 

By annihilating in this equation, a constraint is obtained that determines g^°\ Just as 
in the tokamak calculation, g^ 3 ' = is a solution in the trapped part of phase space. For 
the passing part of phase space, the annihilation operation is {(B/v\\)(-)). Recalling ([IE]) , 
the constraint equation becomes 

= -V 7 r d X I 3 dB cL(°)-A^l. (46) 

V'^'Jo X Jb v\\B-VB Y dt/j J K 1 

We next argue that the S term in A (in (142 p ) may be dropped in ( I46p . The collision 
operator may be written as derivatives and integrals involving v and £ = er\/l — \B, where 
a = sgn(t>ii), so C does not introduce any dependence on x or 7- Thus, the contribution 
to ( 146]) from the dh' /dx term in A is proportional to the leftmost expression in ( 132]) . so it 
vanishes. The constraint equation therefore reduces to 

l B rl \ m (MG + NI\ gigjo" 



°^W C K v -K-^rw)iW (47) 

This equation is identical to the one solved for a tokamak, aside from the constant factor in 
parentheses. 

We now show that the parallel flow associated with the distribution function j\ is consis- 
tent with the form (T33 ]) -( l34l) found in the previous section from a fluid approach. Forming 
/ d 3 vv\\f = riVu using ( I44p and taking g pa g(°\ 



11 B {dip Zen dip 



/MG + NI\ 
V N-tM J K J 



N-iMJo Jo d X (B')Vl - AS' J ' 
where X = J d 3 vv\\g. Noting 

d 3 vQ = irBj2° dX dv—t (49) 



o 



for any Q, the upper limit of the A integral can be changed to 1/B in the integral for X 
because g = in the trapped region. Thus, X varies with position only through the factor 
B, so X has the form of the (VuB) B / (B 2 ) term in (133]) . Next, the A integral in (|4"8]) can be 
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evaluated by moving it inside the B' integral. In this exchange, the range of the A integration 
becomes (0, 1/B') and the range of the B' integral becomes (B, B). When the A integral 
is evaluated, the result is identical to the W term in (I34|) . Consequently, the parallel flow 
(|48p evaluated from kinetic theory has precisely the spatial dependence calculated from fluid 
theory in §M) and (13I]) . 

For the ions in a pure plasma, gi and X can be calculated explicitly using the momentum- 
conserving model collision operator Ci = vC j/ji — /jo m i Ui, ||/^}) where 

*-%k< < 5 °> 

is the Lorentz pitch-angle scattering operator, 

u= ( [ d 3 vf t0 ^u) 1 / <Pvf a w\u (51) 



(52) 



37- 

2nZ 4 e 4 rii In A [erf (x) - # {x)\ 



\l/ (x) = [erf (x) — x (derf (x) /dx)} / (2x 2 ) , erf (x) = (2/s/7r) J x exp (— t 2 ) dt is the error 



function, and x = vj \j2Ti/mi. This model operator captures the dominant effect of colli- 
sions for large aspect ratio. The solution of ( )47|) is performed exactly as for the tokamak 
calculation^, giving 

(MG + NI) m iC dTi (m lV 2 \ av [V* d\' 



where H = H(B 1 — A) is a Heavyside function which is 1 for passing particles and for 
trapped particles. Thus, the ion parallel flow for low collisionality is 

{MG + NI)cB dTj c {MG + NI + W) (d% 1 d Pi \ 
i[l ~ ' h (N - *M) Ze (B 2 ) dip B (N - *M) I # Zen { dip ) ' [ ' 



* 3 /n2\ f l/ * X d\ 

' 1 = iWi TTra (55) 



where 



is the effective fraction of circulating particles. 

When the average parallel ion flow (Vi\\B) is evaluated from (]54p . the W term - the only 
term that reflects the departure from quasisymmetry - does not contribute due to ( 1321) . The 
resulting expression for (VmB) in an omnigenous plasma is identical to the corresponding 
formula for a quasisymmetric stellarator-^. It can be seen that the bootstrap current (j\\B) 
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in an omnigenous stellarator must also be given by the same expression as in a quasisym- 

metric stellarator. To understand this result, first recall that the perturbed distribution 

function of each species is given by f x = g(°' + Ze&if /T + Adf /dip, with A given by (l4"2"j) . 

and the solution of (I47p . The S term in A does not contribute to (V\\B) for the species 

as we have seen due to f[52~j) . Also, g^ must be the same as in a quasisymmetric stellarator, 

since in the equation that determines it, ( H7|) . the departure from quasisymmetry does not 

appear. Thus, (V\\B) for each species is the same as in a quasisymmetric device, and so 

(j\\B) is the same as well. The result is 

. v f t ca (MG + NI) (dpi dp e 2.07Z + 0.88 dT e ^n e dT t \ 
J\\B) = — r — v —— — — — + — n e — - 1.17- 
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can be employed, 
with two Laguerre 



f2 + Z) i N ~ * M ) \dip dip a e dip ' Z dip J ' 

(56) 

where f t — 1 — f c is the effective trapped fraction, a = Z 2 + 2.21Z + 0.75, and Z is 
the ion charge. To obtain ( 1561) . the method of page 207 of Ref. 
using the approximate Spitzer function from appendix B of Ref. 
polynomials. The tokamak expressions for (VfiJB) and (j\\B) can be recovered from the 
omnigenous/quasisymmetric stellarator expressions by simply setting N = 0. 

The expression ( 156]) can also be derived by solving the differential equations for the boot- 
strap current in a general stellarator, derived using the Shaing-Callen moment approach^'^'i 5 -. 
In this approach, the bootstrap current in a general stellarator for a pure Z = 1 plasma is 
given by 

m^H^i^- -^- 1 - 17 "*) (57) 

where the "geometric factor" is 

<G.>4( te) -^>f^M4 (58) 



gi = y/1 — XB, g 2 is the solution of 

B ■ V(g 2 /B 2 ) = BxVip- V(l/B 2 ) (59) 

such that #2 = when B = B, and g^ is the solution of 

B ■ V(g 4 / 9l ) = BxVip- V(1M) (60) 

for any A in the range [0, 1/B] such that g± = when B = B. Equation (|59|) for g 2 closely 
resembles equation (l26|) for U which we solved earlier. Using the solution ( 13T1) . then 

(61) 



1 

92 



N-iM 



(MG + NI) (*--i\- W 
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Equation (!60|) for g± may be solved in the same manner used to find U in ( l26]) -( l3~Tj) . yielding 

VI - \B 1 A j- — rB dB' dti 



MG + NI 
9i ~ N — 7m 



1 + [ 



(62) 



Evaluating (EE}, the W term in #2 and the dh'/dx term in g<± vanish upon flux surface 
averaging due to (j32"]l . Then evaluating the A integral, we obtain (Gb s ) = (MG+NI)/(t,M — 
N), which is precisely the result for a quasisymmetric stellarator. Then ( )57l) for / ( < 1 
reduces to the Z = 1 limit of ( 1561) . proving the two approaches are consistent. 

Equation (15 6j) gives the current density on one flux surface in terms of 7 and G, which 
represent the toia/ poloidal and toroidal current through a suitable surface (times 2/c). By 
writing I and G as the appropriate integrals of the current density, as shown in Appendix 
[C| two ordinary differential equations (1C3|) and ( 1C4[) can be derived that give I and G in 
terms of (j\\B). These two equations, together with ( l56l) . constitute a coupled system of 
equations for the self-consistent current profile. 

The case M = (generalized quasi-poloidal symmetry) is noteworthy, for then the sub- 
stitution of ( 156]) into (IC3j) gives dl/dip = (. . .)/. As the boundary condition for / is that 
it vanishes on the magnetic axis, then the self-consistent profile of is / = 0, and so 
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that the bootstrap current 



(j\\B) = everywhere. Thus, we recover the result of Refs. 
in an M = omnigenous device vanishes. (If Ohmic or RF-driven current is present, there 
will be additional contributions to (j\\B) besides ( 1561) . providing an inhomogeneous term in 
the differential equation for so / and the bootstrap current would become nonzero.) In 
contrast, for any omnigenous plasma with M^O, the contribution from G to the bootstrap 
current is still present. As G has a nonzero boundary condition at the plasma edge, then G 
will generally be nonzero, and so the bootstrap current will also be nonzero. 



IV. RADIAL ELECTRIC FIELD 

In a non-quasisymmetric stellarator, there is usually only one special value of the radial 
electric field E r at each radius for which the electron and ion particle fluxes will be equal. In 
equilibrium, E r must therefore take on this value. (If the electron temperature T e is much 
higher than the ion temperature Tj, a second "electron root" solution may also be possible, 
but we will assume T e ~ Tj, excluding this possibility.) In this section we will determine E r 
for the case of ions in the long-mean-free-path regime. 
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We begin by finding the radial particle flux of each species, allowing general collisionality 
for the moment. By applying (B~ 2 B x Vip ■ (. . .)) to the fluid momentum equation with 
a diagonal pressure tensor, the particle flux is found to be (r • V^) = T m + T E + T d , 
where T m = (J d 3 v fv m ■ Vip), = {ncB~ 2 B x V$ • V^) and T d = c(B~ 2 B x Vip • 
/ d 3 vvC{f}) /(Ze) is the classical flux, which we henceforth neglect. The earlier definition 
(jlOj) with (H2D may be substituted into T m . Using (j4"5j) . the fact that (B ■ VQ) = for any 
single-valued Q, and the fact that A = at A = 1/5, we derive T m = — \J d 3 v Av\\b ■ V/ ) . 



We then substitute in the drift-kinetic equation (1391) . The resulting dfo/dx/j term vanishes 
due to a parity, and the $i term can be shown to cancel the earlier Te term in the total flux; 
this can be done by again applying (jlDjl and noting (B x V^- V($i/-B 2 )) = (V- (§iB~ 2 B x 
V?/>)) = 0. We thereby obtain 



This result holds for any collisionality regime. 

Next, we flux-surface- average the quasineutrality equation V • j — to obtain = 
(j ■ V'lp) = J2 a Z a e(r a ■ Vif>), where a is the particle species. We consider a pure plasma 
with ion charge Z and comparable electron and ion temperatures. In this case, the ion 



contribution to the radial current is larger than the electron contribution by ~ ymj/m e , 
so the leading-order ambipolarity constraint is (Tj • V^) = 0. We henceforth assume all 
quantities refer to ions and drop the subscripts. 

Specializing to the long-mean-free-path regime, we may apply the model collision oper- 
ator from (l50l) - (!52l) . appropriate for large aspect ratio. Using the momentum conservation 
property / d 3 vv\\C{f} = 0, 



We now exploit the fact that {Q dh/dx) = for any Q that is independent of branch and x- 
This fact follows from (fT6|) and fl32l) . Thus, only the terms in (l64"j) that are quadratic rather 
than linear in dh/dx will contribute. For instance, g will not contribute since it is constant 
on a flux surface. Keeping only the nonvanishing terms, we may write (r • V^) = r\ + T 2 
where 




(63) 





(64) 




(65) 
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/ 



ihSvC \;^ )h (66) 



and 

^=-(/^^)"7 d3 "I wi|S (67) 

is the part of u that depends on Applying ( 1491) . it can be seen that I\ and w x are both 
proportional to the integral oc J °° dvv^vdfo/ dip. Therefore, T2 and the total radial flux are 
proportional to the same factor. The ambipolarity condition can thus be written 

^ dh'dh"\ , 4 df 



°=< r ' v « K U *<**jy. ivv "t (68) 

where the single and double primes refer to the integration variables in the S factors. In a 
quasisymmetric or axisymmetric plasma, dh/dx = and so the equation is automatically 
satisfied. However, in a non-quasisymmetric omnigenous field, the x integral in fl68l) is 
generally nonzero, so the v integral following it must vanish. This condition may be written 



f°° A - 2 

= / dx x e x v 
Jo 



1 dp Zed$ ( 2 5\ 1 dT 
pd^ + ~T^hf + V ~2)T~d4 



(69) 



which may be solved for the radial electric field d^o/dip. Reinstating the species subscripts 
for completeness, the result is 

d% Ti ( 1 dpi 1.17 dTA 



dip Ze \ pi dip Ti dip 

where 



(70) 



1.17= --(^y dxx 4 e- x uj J dxx 6 e~ x v (71) 

is the same numerical constant that appears in the banana-regime tokamak ion flow and 
in ( 15^1) . This electric field differs from previously known results for a non-omnigenous stel- 
larator. For example, if the main ions in a non-omnigenous device are in the \ jv regime of 
collisionality, the relation f lTD]) holds but with 1.17 replaced by 2.37 (see e.g. equation (35) 



of Ref. |25j). 

Assuming the temperature scale length is not dramatically shorter than the density scale 
length, ( 1701) gives an inward electric field. Physically, this field develops to electrostatically 
confine the ions, reducing their radial flux to the much smaller level of the electron flux. 

As the departure from omnigenity increases, (1701) will become inapplicable before the 
formulae for the flow and current do, since the particles with nonzero average radial drift 
will contribute strongly to the radial particle flux but not to the parallel transport. 
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V. CONSTRUCTING STELLARATOR-SYMMETRIC OMNIGENOUS 
FIELDS 



In this section we will give a construction for B(9, £) field strength patterns that both 
satisfy all the omnigenity conditions and also satisfy stellarator symmetry^, meaning the 
invariance of B under the replacements 9 — > —9, £ — > —( . A construction of omnigenous 
B(9, was given previously in Ref. Q, but that procedure generally produces fields without 
stellarator symmetry, whereas every stellarator experiment to our knowledge does possess 
stellarator symmetry (aside from small error fields). 

The fields we will construct are designed to have the following omnigenity properties: 1) 
The B contours will all link the torus, 2) the B contour will be straight, and 3) dA^/dx = 
and dAg/dx = 0. Here, A^ and A# are the same quantities discussed following (fT4|) : the 
separations in £ and 9 between the pair of points on opposite branches of a field line but at 
the same B. The other omnigenity properties from section [III will then follow automatically. 
For example, by applying d 2 /dBdx to ( |T4l) . then ( |T2l) follows. 

We will now present the construction, and we will verify afterward that it indeed produces 
fields that are omnigenous and stellarator- symmetric. The construction is given in terms of 
new angles 9 and (, such that consecutive maxima of B lie on the constant-^ curves ( = 
and ( = 2tt. We also employ an effective rotational transform 1, equal to d9/d( along the 
field. To construct an omnigenous field with nonzero iV and with N P N toroidal periods, 
the new quantities are related to the original quantities by 9 — 9, £ = (N( — M9)N P , 
and I = i/[{N — iM)N p ]. To construct an omnigenous field with (M,N) = (1,0) and N p 
toroidal periods, then instead 9 = N p (, ( = 9, and I = N p /t. In either case, a stellarator- 
symmetric B is one invariant under 9 — > —9 and ( — > — £. Also, dA^/dx = is equivalent 
to dA^/dx = 0, where is defined just as for A^ but using ( in place of £. 

There are several inputs to the construction. First, we may pick B and B. Second, we 
choose a function D(x) which will turn out to be closely related to A^(B). The function 
D(x) must be defined on the domain [0, n] and must satisfy -D(O) = it and D(tt) = 0. Lastly, 
we choose a function s(x,y) which will represent the (^-variation of the B contours. We 
require s(x, y) to be both odd in y and 27r-periodic y (i.e. the Fourier series for s contains 
sin(n?/) terms but no cos(ny) terms or y-independent term.) The input x ranges over [0, 7r], 
and we require s(0, y) = for all y. 
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FIG. 8. (Color online) Contours of rj(6, () for the field of Figure [3] One branch (77 < tt) is shaded 
while the other branch (77 > tt) is unshaded. Field lines are parallel to the arrow. P_ and P+ are 
defined preceding ([75]). 

It is then useful to introduce a new coordinate rj which resembles B but which is different 
in the two branches. Specifically, we define 77 G [0, 2tt] through the relation 



B = B(l + e + e cos 77) 



(72) 



where e = (B - B)/(2B). Notice from Q7JJ) that 7/ = along ( = (where B = B) and 
77 = 27r along £ = 2tt (where 5 again rises to B). While B contours coincide with 77 contours 
in the (#,C) plane, 77 lies in the range < 77 < n on one branch, while 77 lies in the range 
< r] < 2tt on the other branch. 
We then compute ( as follows: 

7r - s (77, 6 + ZD(rj)) - D(rj) if < 77 < tt 



CM) = < ~ x ( 73 ) 

7T + s (27T -7], -6 + lD{2h - 77)) + ,D(27r - 77) if 7r < 77 < 2vr. 



Numerical root-finding is next used to compute the inverse map r](6,(), an d finally B is 
calculated by ( 172]) . 

Figures E] and S] show omnigenous stellarator-symmetric fields constructed using the above 
procedure. For Figure [31 the parameters are chosen to resemble HSX-: M — 1, N — 4, 
iVp = 1, * = 1.05, and e = 0.072. The numerical functions used are s(x,y) = 0Axsin(y) and 
D(x) = it — x. Figure [8] shows the associated r](9, () function. For FigureHJ M = 1, N = 0, 
N p — 3, * = 1.62, e = 0.1, s(x, y) = 0.15xsin(7/) and D{x) = n — x. 
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We now prove that the B resulting from the above construction is stellarator-symmetric 
and omnigenous, beginning with stellarator symmetry. Thinking of B as an independent 
variable in place of (, we can restate the definition of this symmetry as follows: when 9 — > —9, 
and B remains constant, and the branch is reversed, then ( must go to — £. Reversing the 
branch at constant B is equivalent to the replacement 77 — ^ 27r — 77. Adding factors of 2w 
to keep 9 and ( within the range [0, 2tt], we can therefore write the stellarator symmetry 
criterion as 

CM) = 27r-C(27r-77, 2tt-0). (74) 

If 77 < 7r, then the left-hand side of this equation is given by the top line of (1751) . and the 
right-hand side of (1741) is given by the bottom line of (173]) . It can be immediately verified 
that ( 1741) is indeed satisfied. Similarly, if 77 > it, then the left-hand side of ( 1741) is given by 
the bottom line of (1731) . while the right-hand side of (1741) is given by the top line of (1731) . 
and the satisfaction of (174)) is immediate. 

It remains to show that A^ is independent of field line. To see this, consider two points P + 
and P- on the same field line and at the same B but on opposite sides of B. Figure |8] shows 
two such points. Let ( + , 7? + ) describe the point P + with r/ > tt and let (§_, r/_) 
describe the point P_ with rj < 7r. Applying (1731) to P_, 

C- = tt - 8 (V, L + W (r/„)) - D fo_) . (75) 

Applying (1751) to P + , noting 77+ = 2n — r\_ and 9 + = 9_ + — C-): an d recalling s is odd 
in its second input, then 

C+ = tt - s (r/_, 0_ + t(C + - C-) - ^(r/„)) + D(77_). (76) 

Comparing ( 1751) and ( 1761) . it can be seen that a consistent solution of the system is ( + — = 
2D{j]_). Therefore = ( + — is independent of 9 for each B, so the magnetic field is 
omnigenous. 

Given that * is rarely much larger than 1 in experiments, it is difficult to construct = 
omnigenous fields that depart strongly from quasisymmetry. This is because the B contours 
in a N = quasisymmetric field are already nearly parallel to the field lines when t/N p < 1. 
Even a slight curvature of the B contours would result in points of tangency to the field 
lines, and as discussed in Section III B] the B contours and field lines can never be tangent 
in omnigenous fields. 
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VI. CONCLUSIONS 



The limit of a single-helicity (i.e. quasisymmetric) field is a useful point of reference for 
insight into stellarator physics, for in these effectively 2D fields, analytic calculations can 
be carried out more completely than in a general stellarator. In the preceding sections, we 
have shown that omnigenous plasmas ought to be another such point of reference. Many 
geometrical and physical properties of omnigenous plasmas are either identical to the associ- 
ated quantity in a quasisymmetric plasma, or else obtained by the addition of one term that 
is strongly constrained. While every quasisymmetric field is omnigenous, in practice not 
every omnigenous field is quasisymmetric. Therefore a broader attention to the larger class 
of omnigenous fields may allow better optimization for other criteria such as smaller aspect 
ratio or coil simplicity. Furthermore, it is omnigenity and not quasisymmetry that is the 
necessary condition for confinement of alpha particles. Any viable fusion reactor will need 
to be nearly omnigenous in order to prevent damage to the first wall from unconfined al- 
phas. Indeed, alpha particle confinement time is routinely used as an optimization criterion 
in stellarator design codes. As alpha confinement becomes increasingly important in future 
reactor-relevant experiments, stellarator desi gns can be expected to more closely approach 
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gives one example of a nearly omnigenous 



omnigenity. The recent design study in Ref. 
device. 

In Section [III novel proofs were given of various geometric properties of omnigenous 
magnetic fields. The B contours must link the flux surface toroidally, poloidally, or both, 
so no isolated extrema of B on the flux surface are allowed. The field lines can never 
be tangent to the B contours. The curve of B (the maximum of B) must be straight in 
Boozer coordinates. The variation of B on the two branches is constrained by (Tl2|) . The 
integers M and N usually defined by the quasisymmetry relation B = BiMQ — NQ may 
be generalized to any omnigenous field by redefining them as the number of times each B 
contour encircles the plasma toroidally and poloidally. In ( TTTT) . it was shown that the ratio 
(B x V-B ■ V'ip)/B ■ V-B, which arises repeatedly in transport calculations, may be written 
as a sum of a quasisymmetric part and a departure from quasisymmetry, dh/dx- 

This same division into quasisymmetric and non-quasisymmetric components also arises 
in physical quantities. The Pfirsch-Schliiter flow and current ( |33|) and ( |36|) are given by 
the same formulae as in a quasisymmetric device, with the addition of the W term, which 
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is an integral of B along field lines. In contrast, the Pfirsch-Schliiter flow and current in 
a general stellarator can only be expressed in terms of the solution of a partial differential 
equation like ( I26p . The distribution function in a low-collisionality omnigenous plasma is 
given by the distribution function in a quasisymmetric plasma of the same helicity, plus 
the term —S dfo/dip. This term does not contribute to the average parallel flow (VmB) 
and bootstrap current (j\\B), so these quantities are given by exactly the same expressions 
as in a quasisymmetric plasma. For each of the quantities discussed above, formulae for a 
quasisymmetric plasma can be recovered by setting dh/dx — > 0, and tokamak formulae can 
be recovered by also setting iV — > and M — > 1. 

In general, a self-consistent current profile may be found by calculating the current den- 
sity in term of the total current using kinetic theory, and then writing the current density 
as the appropriate derivative (IC3j) -f lC4|) of the total current. When M — 0, correspond- 
ing to generalized quasi-poloidal symmetry, the resulting self-consistent configuration has 
no bootstrap current. Such configurations have several desirable properties owing to the 
minimization of current driven by pressure. The magnetic field would remain omnigenous 
over a larger range of plasma pressure. Also, the equilibrium (3 limit may be less severe and 
current- driven instability may be reduced. 

The property of omnigenous plasmas that is most fundamentally different from quasisym- 
metric plasmas is the radial electric field. Omnigenous plasmas that are not quasisymmetric 
are not intrinsically ambipolar. The radial electric field E r is therefore determined by low- 
order neoclassical physics, so it is possible to solve for E r using the ambipolarity condition. 
The closed-form expression f JTOj) gives the resulting radial electric field for low collisionality. 
The formula is independent of the details of the magnetic field geometry. 

In every stellarator experiment to date, detailed numerical neoclassical calculations^ find 
a, 1/u regime of radial transport to exist, indicating that departures from omnigenity are 
non-negligible for radial transport in realistic experiments. Thus, the results in Section IIVI 
are likely to be applicable only for benchmarking numerical codes, since in codes the B(6, Q 
pattern can be made more perfectly omnigenous than in a true 3D equilibrium. However, 
unlike the radial transport, the parallel flow and current are carried by the bulk of the 
distribution function, not by a small fraction of trapped particles. Therefore, the formulae 
derived in this paper for the flow and current are robust to small departures from omnigenity. 
This resilience is aided by the fact that any 1/u part of the distribution function would be 
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even in sgn(wi|) and so it would carry no flow. 

In conclusion, omnigenity is a useful ideal limit for gaining insight and benchmarking 
codes: omnigenity is experimentally relevant, it is more inclusive than quasisymmetry, and 
yet it allows explicit analytic results to be derived that are nearly as concise as results for 
quasisymmetric or axisymmetric plasmas. 



ACKNOWLEDGMENTS 

This work was supported by the US Department of Energy through grant DE-FG02- 
91ER-54109 and through the Fusion Energy Postdoctoral Research Program administered 
by the Oak Ridge Institute for Science and Education. 



Appendix A: Additional Proofs 

In this appendix we will first prove the equivalence of the two definitions of omnigenity: 
1) the bounce-averaged radial drift vanishes, and 2) the longitudinal adiabatic invariant is 
a flux function. Then, we will derive (I12p . 

Begin by considering a general (not necessarily omnigenous) stellarator equilibrium. We 
then write Vm-Vip = (v||/Q)Vx (v\\b)-'Vt/j = (w||/f2)V- (v\\/B)B x , where the gradients 
are evaluated at fixed A and v. The formula for the divergence in a general coordinate system 
is then applied, using the coordinates (ip,a,(), where a = 9 — is a field line label, and 
9 and ( are any straight-field-line poloidal and toroidal angles (so B = x Va). Noting 
that the inverse Jacobian is 1/y/g = B ■ VC and that B x \7ip ■ VC = —I/yfg, we thereby 
obtain 

rnr I f) \ I iin \ I f) \ I Tihi \ 

(Al) 



v, 



TTtC 

V^ = —v {l b-VC 



d_ 

da 



b ■ VC 



d_ 



Iv\\ 
~B 



Next, we define the bounce average, which for any quantity A is A = r 1 Yla ° Ic-( v \\^ ' 
VQ^AdC, where r = 2 J c C _ + (|t'|||6- VC)" 1 ^, 

a = sgn(ui|), and and (+ are the two bounce 
points (at which i>m =0). The integrals in the bounce average are performed at constant a. 
Applying the bounce average to flAip gives v m 

Even though and ( + depend on a, it is valid to pull the dj da derivative in front of the 
integral in this result because the integrand vanishes at these endpoints. Now, consider the 
longitudinal invariant for trapped particles: J (ip, a, v , A) = $vud£, where the integration 
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is again carried out along a full bounce. Using d£ = d(/b ■ V( in this definition, then we 
obtain 

mc dJ 



*"'' V ^I^a (A2) 
Consequently, the bounce-averaged radial drift vanishes if and only if the longitudinal in- 
variant J is a flux function, and so the two definitions of omnigenity are equivalent. 

Now we move on to the proof that (d/dx) S 7 l/b ■ V£> = in an omnigenous field^. We 
first note that J = 2vJ2jl fjj X {b ' VB) _ Vl — XBdB. Applying d/dx an d requiring J to 
be constant on a flux surface gives 

0= I" dBy/x^B S (B) (A3) 

J B 

where x = 1/A and S (B) = (d/dx) J2yl/(b • VB). Equation (1A3j) is true for any x in the 
interval (B, B). We now divide AA3J) by y/y — x, where y is any value in the same interval, 
and we integrate over x from B to y. Interchanging the order of the x and B integration, the x 



integral becomes /jj dx^J (x — B)/(y — x) = (vr/2) (y — B), giving = J| dB (y — B)S (B). 
Differentiating twice with respect to y gives = S (y). We had allowed y to be any value in 
the interval (B, B), so S must vanish everywhere, proving the theorem. 

Appendix B: Quasisymmetry in various coordinate systems 

In this appendix, we prove that B has a single helicity in Boozer coordinates if and only 
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if B has a single helicity in Hamada coordinates. While one proof can be found in Ref. 
here we prove a more general result, that symmetry is equivalent for any straight-field-line 
coordinate system in which the Jacobian is proportional to some power of B. 

A number of identities must be demonstrated before the main theorem is proved. Begin 
by considering two sets of coordinates, (6 X , ( x ) and (9 y , Q), which are not necessarily Boozer 
or Hamada coordinates. If the the coordinates are straight-field-line coordinates, then B = 
Vip x V9 X + *VCc x = Vifj x V6 y + *VC y x where 2nip is the toroidal flux. Suppose 
the transformation from one system to the other is written as 9 y — 6 X + T (ip, 9 X , ( x ), ( y = 
(x + Q (ip, 9 X , Cx) where J 7 and Q are periodic in both the poloidal and toroidal angles. Then 

V^xVJ + *V0 x W = 0. (Bl) 

The V9 X component of this equation tells us dF/d(, x = tdQ/d( x , so upon integrating, 
T — i>Q + y (ip, 9 X ). Here and throughout this appendix, d/d( x holds 9 X fixed, d/d9 x holds 
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C, x fixed, d/d( y holds 9 y fixed, and d/d9 y holds Q fixed. The V( x component of (1B1|) implies 
dJ-'/dOx = idQ/d9 x , so J 7 = tQ + w(ip,(x)- By comparing these two relations between J 7 
and Q, T must equal *£/ plus a flux function, so 

9 y = 6 x + *g + A(rl>) and C y = C* + Q (B2) 

for some flux function *4 (V0- 

Now suppose the Jacobian for the (6 x ,Cx) coordinates is proportional to B x , that is, 
VV' • V# x x VCx = {ip) B x for some flux function A x (ip). Hamada coordinates have x = 
and Boozer coordinates have x — 2. The flux surface average of a quantity Q is 

1 r 2n r2n Q 

K ' JO JO ■ Vt^ X VCa; 

where V (-0) is the volume enclosed by a flux surface, and the prime denotes d/dip. Observe 
that (B x ) =4tt 2 /(V'A x ), so 

Ait 2 B x 

^■V9 x xV( x = —j^- y (B4) 

Now suppose the two straight-field-line coordinate systems (9 X , ( x ) and (9 y , ( y ) have Ja- 
cobians proportional to B x and B y respectively. Then from ( 1B4I) . 

V0 ■ x VC, = ^fr^yW ■ V9 X x VCr- (B5) 



Next, we apply the chain rule to Q from (IB2j) : 

dg_^d9,dG_ dQdG_ 

89 y 89 y d9 x 89 y d( x [ ' 

By applying 0/00„ to dm}, we find 1 = 89 x /d9 y + idG/d9 y and = d( x /d9 y + dG/d9 y , so 
( 1B6P implies 

89 y { h 89 y ) 89 x d9 y d( x 1 ' 

,dG dg\dg_ 8Q 

{ 1 + dZ + i d9-)dT y -d9- x - m 
Now, apply B ■ V to (jB2]) to obtain x V9 y ■ V( y = x V6» x • VCx + B ■ VQ. Noting 
(!B5l) . then 

og og 
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Rearranging, 



Substituting this expression into (1B8j) then gives 

(B x ) Bv dg 



eg 



(BIO) 



B x (By)d6 v d6 x 

Repeating the steps from flB6j) - flB8j) but applying d/d6 y instead of d/d( y , we can similarly 
derive 

(b x ) By og _ og 
b* (By)~dc y ~~dc~ x 

Next, applying d/d6 x to flBQl) . 



(Bll) 



(y 



x) ^l B y^ 1 — = (— + ^~ 
! (By) de x \d( x dd x ) w x 



(B12) 



Recalling that B ■ V = W x V6 X ■ V(x [(d/d( x ) + * (d/dO x )], then flBT2l) is equivalent to 

dg 4vr 2 (y 



B - V 



(B13) 



where we have also applied (1B4[) . We could have applied d/d( x to (1B9[) instead of d/d9 x , 
and so it is also true that 

B . v « = 4^fe-x) 9£ 

aCx v (s»> ac. 1 ' 

We are finally prepared to begin the main proof. Suppose B has only a single helicity in 
the (0 x ,(x) coordinates: B = B (M8 X — N( x ) for some integers M and N, or equivalently, 
NdB/d6 x + MdB/dC x = 0. Then from flBT3|) - (lBT4l) . B ■ V (Ndg/d6 x +MdQ/d( x ) = 0. 
It follows that N dg I d6 x + M dg / d( x = S (ip) for some flux function S (-0) . Integrating this 
result in 9 X and ( x from to 27r in both variables, we find S must be zero. Then applying 
dBTOl and (IBTTjl . 



Finally, we form 



N 



dB 



M 



dB 



N 



N 



N dg/d9y + M dg/dCy = 
f d6 x dB d( x dB" 



M v d6 x 



og_ 



09 y d(a 

dB 
df x ~ 



-M\ 

dg dB' 
d9 y dC x , 



d9 x dB d( x dB 



d( y d6 x 



M -i 



d(y d( x 

dg dB 

d( y d6 x 



(B15) 



(B16) 





dB\ 


. 9(y_ 


dx~ x ) 



The first equality above is the chain rule, and to get the second line we have used the d/d6 y 
and d/d( y derivatives of ( 1B2[) . The last line of ( 1B16[) vanishes due to (1B1 5|) . and so 

(B17) 



^ T dB t 8B 
N— + M— 

d0 x d( x 







dB flB 
N— + M- 



0. 



d9y ' " d(y 

The right equality in (IB 171) also implies the left one, since x and y were arbitrary in the 
proof. For the specific case of x = and y = 2, then B has a single helicity in Boozer 
coordinates if and only if B has a single helicity in Hamada coordinates. 
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Appendix C: Current in a general stellarator 



Here we calculate several relations which are satisfied by the current in any MHD equilib- 
rium with nested toroidal flux surfaces. We begin by noting that the perpendicular current 
is j ± = c (dp^/dip) B~ 2 B x Wip, where is the sum of the pressures of each species and a 
flux function. 

Next, recall that the coefficient / (ip) in the covariant Boozer representation ([3]) equals 
2/c times the toroidal current inside a flux surface. Therefore / (ip) = (2/c) / d 2 a ■ j, where 
the surface integral is performed over a constant-^ cross-section of the plasma, a surface 
which covers the region from magnetic axis out to the flux surface ip. The area element is 
d 2 a = dipdO (Vip ■ V0 x VC)" 1 VC Using j = (j\\/B)B + j ± , Q2J, and ©, we obtain 

2 r "#'f- c /^r^+r^v (CD 



c Jo \ dtp Jo B 2 Jo B 

where everything in parentheses is evaluated at ip' rather than ip. We next integrate f ICTI) 
over all (. Recall that the flux surface average in Boozer coordinates can be written 



S^de^dai/B' 



and observe that (B 2 ) = An 2 / d6 J 27r d( (1/B 2 ). Then differentiating in ip, we can 
write 

dJ_ 47r/dps = 4tt (j\\B) 

dip (B 2 ) d^P c (B 2 ) ' 1 j 

The boundary condition for this equation is / = at the magnetic axis. 

An analogous ODE for G(ip) can be derived by repeating the above analysis with a 
constant-6 1 surface: 

dG AirGdpy 47r *(j||£) 

dip (B 2 ) dip c (B 2 ) ' 1 } 

The boundary condition for G is that it must go to its vacuum value at the plasma edge. 

If another equation for (j\\B) in terms of I and G can be obtained from kinetic theory 
(as we have done for an omnigenous stellarator in (|56l) ). then this equation can be used with 



f)C3j) and (IC4I) to calculate a self-consistent current profile. 
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